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Abstract. The Square Kilometre Array (SKA) is a planned multi purpose radio telescope with a collecting area 
approaching 1 million square metres. One of the key science objectives of the SKA is to provide exquisite strong- 
field tests of gravitational physics by finding and timing pulsars in extreme binary systems such as a pulsar-black 
hole binary. To find out how three preliminary SKA configurations will affect a pulsar survey, we have simulated 
SKA pulsar surveys for each configuration. We estimate that the total number of normal pulsars the SKA will 
detect, using only the 1-km core and 30 minutes integration time, is around 14 000 normal pulsar and 6 000 
millisecond pulsars. We describe a simple strategy for follow-up timing observations and find that, depending on 
the configuration, it would take 1-6 days to obtain a single timing point for 14 000 pulsars. Obtaining a single 
timing point for the high-precision timing projects of the SKA, will take less than 14 hours, 2 days, or 3 days, 
depending on the configuration. The presence of aperture arrays will be of great benefit here. We also study the 
computational requirements for beam forming and data analysis for a pulsar survey. Beam forming of the full field 
of view of the single-pixel feed 15-m dishes using the 1-km core of the SKA requires about 2.2x10^^ operations 
per second. The corresponding data rate from such a pulsar survey is about 4.7x10^^ bytes per second. The 
required computational power for a deep real time analysis is estimated to be 1.2x10^^ operations per second. 
For an aperture array or dishes equipped with phased array feeds, the survey can be performed faster, but the 
computational requirements and data rates will go up. 
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1. Introduction 

There are, arguably, no other astronomical objects whose 
discovery and subsequent studies provides more insight 
in such a rich variety of physics and astrophysics than 
radio pulsars. Pulsars fi nd their applications (see e.g. 
Lorimer fc Kramen 12005! ) in the study of the Milky Way, 
globular clusters, the evolution and collapse of massive 
stars, the formation and evolution of binary systems, the 
properties of super-dense matter, extreme plasma physics, 
tests of theories of gravity, the detection of gravitational 
waves, and astrometry, to name only a few areas. Indeed, 
many of these areas of fundamental physics can be best - 
or even only - studied, using radio observations of pulsars. 

To harvest the copious amount of information and sci- 
ence accessible with pulsars, two different types of obser- 
vations are required. Firstly, suitable pulsars need to be 
discovered via radio surveys that sample the sky with high 
time and frequency resolution. Depending on the partic- 



ular region of sky to be covered (e.g. along the Galactic 
plane vs. higher Galactic latitudes), different technical re- 
quirements may be needed. Secondly, after the discov- 
ery, a much larger amount of observing is required to 
extract most of the science using pulsar timing obser- 
vations, i.e. the regular high-precision monitoring of the 
pulse times-of-arrival and the pulse properties. Again, de- 
pending on the type of sources to be monitored (e.g. young 
pulsars in supernova remnants vs. millisecond pulsars in 
Globular clusters) the requirements are different. 

In the past 40 years, astronomers have impressively 
demonstrated t he potential of pulsar s earch and timing ob- 
servations (e.g. Hulse fc Tavlod l 1975V However, the next 
10-15 years promise to revolutionise pulsar astrophysics in 
a way that is without parallel in the history of pulsars or 
radio astronomy in general. This revolution will be pro- 
vided by the next generation radio telesc ope known as the 
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"Squ are-Kilometre- Array" (SKA) (e.g. 
20061 ). The SKA will be the largest telescope ever built, 
with a maximum baseline of 3 000-f km and about a factor 
of 10-100 more powerful (both in sensitivity and survey 
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speed) than any other radio telescope. The SKA will be 
particularly useful for pulsars and their applications in as- 
trophysics and fundamental physics. Due to the large sen- 
sitivity of the SKA, not only a Galactic census of pulsars 
can be performed, but the discovered pulsars can also be 
timed more precisely than before. As a resu lt, the pulsar 
scienc e ena bled by the SKA and described bv lCordes et al 
(|2004[ ) and lKramer et al.l (|2004l ) includes: 



A 3 000 15-m dishes with a single-pixel feed, a sensitivity 
of 0.6 SKA units, Tsys=30K and 70% efficiency cover- 
ing the frequency range of 500 MHz to 10 GHz. 

B 2 000 15-m dishes with phased array feeds from 
500 MHz to 1.5 GHz, a sensitivity of 0.35 SKA units, 
a field of view (FoV) of 20deg2, T^y^^SSK and 70% 
efficiency and a single-pixel feed from 1.5 to 10 GHz, 
with rsys=30K. 



— Finding rare objects that provide the greatest oppor- 
tunities as physics laboratories. Examples include: 

— Binary pulsars with black hole companions, which 
enable strong-field tests of General Relativity. 

— Pulsars in the Galactic Centre, which will probe 
conditions near a 3 x 10^ Mq black hole. 

— Millisecond pulsars that can act in an ensemble 
(known as Pulsar Timing Array) as detectors of 
low-frequency (nHz) gravitational waves of astro- 
physical and cosmological origin 

— Millisecond pulsars spinning faster than 1.4 ms that 
can constrain the equation-of-state of super-dense 
matter. 

— Pulsars with translational speeds in excess of 
lO'^kms^^ that can probe both core-collapse 
physics and the gravitational potential of our 
Galaxy. 

— Rota ting radio transients (see McLaughlin et al. 3. Pulsar surveys 
l2006l) 



C A combination of a dense aperture array (AA) with a 
FoV of 250 deg^, a sensitivity of 0.5 SKA units, cover- 
ing the frequency range of 500 to 800 MHz and 2 400 
15-m dishes with a single-pixel feed covering the fre- 
quency range of 800 MHz to 10 GHz, a sensitivity of 
0.5 SKA units, rsys=30K and 70% eflfciency. 



We assume that 20% of the elements will be placed 
within a 1km radius and 50% within a 5 km radius. The 
sparse aperture array below 500 MHz will be ignored, as 
we will not consider using it for pulsar surveys or timing in 
this paper. We will refer to the dishes from configurations 
A and B (which have a lower frequency limit of 500 MHz) 
as 500-MHz dishes and we will refer to the dishes from 
configuration C (which have a lower frequency limit of 
800 MHz) as 800-MHz dishes. 



— Intermittent pulsars (see Kramer et al. 2006^ 

— Understanding the advanced stages of stellar evolu- 
tion. 

— Probing the interstellar medium of our Galaxy. 

— Understanding the pulsar emission mechanism. 

The aim of this paper is to investigate, for differ- 
ent configurations of the SKA and different computation 
power and data rates, a variety of different methods to 
achieve the best results for pulsar searches and pulsar tim- 
ing. We approach this problem by using Monte Carlo simu- 
lations based on our current understanding of the Galactic 
pulsar population. In §2, we summarise the various config- 
urations for the SKA. In §3, based on current population 
synthesis results, we investigate the SKA survey yields as 
a function of configuration and sky frequency. In §4, we 
describe a simple strategy for optimising follow-up timing 
observations. Computational requirements for the pulsar 
surveys are discussed in §5, from which we make specific 
recommendations in §6. Finally, in §7, we summarise our 
results. 

2. SKA configurations 

We will express different sensitivities as a fraction of one 
'SKA unit' defined as 2 x 10'^ m^ K~^. Following the pre- 



One of the main science goals for the SKA is to find the 
majority of the pulsar population in the Galaxy. These 
should include binary and millisecond pulsars, as well as 
transient and intermittent sources. A pulsar survey con- 
sists of two parts. The first part is using the telescope 
to observe parts of the sky with a certain dwell time. 
The total amount of telescope time depends on the FoV 
of the telescope, the chosen dwell time of the observa- 
tion and the amount of sky that is searched. The sec- 
ond part is the analysis of the observations, which can be 
performed either real time, or off-line. Depending on the 
amount of data that needs to be analysed and especially 
if acceleration se arches are employed to se arch for binary 
pulsars (see e.g. iLorimer fc Kra^leIll2005^ ■ this can be a 
very computationally expensive task. Real time processing 
has been done in the past, but any major pulsar survey 
has always had an o ff line processing component (see e.g. 
Lorimer et al. 20061). 



liminary specifications given bv lSchilizzi et al.l \2OOlh . we 
assume the SKA to consist of a sparse aperture array of 
tiled dipoles in the frequency range of 70 to 500 MHz and 
above 500 MHz it consists of one of the following three 
implementations: 



3.1. Technical considerations 

The FoV of a survey has a maximum that depends on the 
characteristics of the elements used. In the case of circular 
dishes with a one pixel receiver, the FoV can be approxi- 
mated as (A/Ddish)^ steradians, where A is the wavelength 
at which the survey is undertaken and -Ddish is the dish 
diameter. However, by placing a phased array feed a t the 
focal point of the dish (see e.g. Ilvashina et al.l 120041 ). the 
FoV can be extended by a factor which we denote ^. The 
total FoV then no longer depends on A. In the case of the 
Aperture Array (AA), the FoV of the elements is about 
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half the sky. The actual size of the FoV that can be ob- 
tained will be limited by the available computational re- 
sources. This is because the signals from the elements will 
be combined coherently, resulting in 'pencil beams', the 
size of which scales with where Z?oi is the distance 

between the furthest elements that are used in an obser- 
vation. To obtain a large synthesised FoV it is therefore 
necessary to restrict the beam forming to using only the 
elements in the core of the telescope (so that = i'corc), 
which will, however, reduce the sensitivity. 

3.2. Survey simulations 

To find out how the SKA pulsar survey performance de- 
pends on different SKA configurations, we simulated SKA 
pulsar surveys for different collecting areas and different 
centre frequencies. 

3.2.1. Simulation method 



The s imulations were performed following iLorimer et al 
and using their Monte Carlo s imulation pack- 



ag^. In their study, iLorimer et al. I (l2006l) used the results 
from recent surveys with the Parkes Multibeam system 
to derive an underlying population of 30,000 normal pul- 
sars with an optimal set of probability density functions 
(PDFs) for pulsar period (P), 1400-MHz radio luminos- 
ity (L), Galactocentric radius (i?) and height above the 
Galactic plane (z). We make use of these results in our sim- 
ulati ons described below w hich adopt the PDFs of model 
Lorimer et aP (|200d ). 



Our simulation pr o cedur e begins, following the find- 



ings of Lorimer et al. I (|2006h . by generating 30 000 nor- 
mal pulsars which beam towards the Earth. Each pul- 
sar is assigned a value of P, L, R and z based on 
the assumed PDFs. Their initial positions in the plane 
of the Galax y follow the spiral-arm modeling procedure 



described by iFaucher-Giguere fc Kasp: 
trinsic pulse width of each pulsar 



(|2006h . The in- 

foUows the power 



la w relationship w ith sp in period given in equation (5) 



of ILorimer et al.l (|2006l ) (see also their Figure 4). To 



compute the expected DM and scatter broadening ef- 
fects o n each pulse, we use t he NE2001 electron density 
model ( Cordes fc Lazicill2002 ). Note that, since we are pri- 
marily concerned with the distant population of highly 
dispersed pulsars in these simulations, we do not attempt 
to account for interstellar scintillation. Finally, to allow 
us to extrapolate the 1400-MHz luminosities to other sur- 
vey frequencies in the next step, we make the reasonable 
assumption t hat pulsar spectra c an be approximated as 
a power law (jLorimer et al.lll995[ ) and assign each pulsar 
a spectral index drawn from a normal distribution with 
mean of -1.6 and standard deviation 0.35. The same pro- 



^ The package psrpop can be obtained from 
|http : // psrpop. sour ceforg e .net [ Additional software to 
simulate and analyse a series of surveys can be obtained from: 
|http : //www, jb.man. ac .uk/~rsniits/SKA. tgz| 



cedure was used to obtain a millisecond pulsar popula- 
tion, but with pulsar period and pulsar luminosity dis- 
tributions from Cordes fc Cherno and assuming 
30 000 p otentially detect able millisecond pulsars in the 



Galaxy (jLvne et al.lll998D 



Given the above model for the underlying population, 
which is consistent with the results of current-day surveys, 
we proceed to apply it to different configurations of the 
SKA. For each configuration, we compute the observed 
flux density of each pulsar and compare it with the limit- 
ing flux threshold, which depends on the sky temperature 
(Tsky) and the observed width of the pulsar profile. For 
the AA we reduced the sensitivity by the cosine of the 
zenith angle of the source, as the AA is sensitive only to 
the component of the radiation projected onto the plane 
of the AA. In addition, because the core of the SKA will 
be located near a declination of -30°, the sky visible to 
the SKA was limited to a maximum decli nation of 50°. 

jiooi). 



As discussed by ILorimer et al 



we make the 

important distinction between the observed (i.e. post- 
detection) pulse width which is related to the pulsar's in- 
trinsic pulse width by the quadrature sum of contributions 
from samplin g, dispersion s meari ng and scattering (see 
equation 4 of Lorimer et al.l ( 2006h ). To scale the scatter- 
broadening time to arbitrary SKA surve y frequencies, we 
adop t a frequency power law index of -4 (jCordes fc Lazio 



20021 ). We consider a pulsar to be detected by a given con- 



figuration if its flux density exceeds the threshold value, 
and its observed pulse width is less than the spin period. 

Initially, we determined the optimal centre frequencies 
for a survey in the Galactic plane (defined as |1| < 45° 
and |b| < 5°) and an all-sky survey excluding the Galactic 
plane. This was done by comparing the detected number 
of normal pulsars by simulating surveys with sensitivi- 
ties in the range 0.1 to 0.5 SKA units (in steps of 0.1 
SKA units) and with centre frequencies in the range 0.4 
to 2 GHz (in steps of 0.25 GHz). For each sensitivity and 
for both survey types, the frequency that provided the 
largest number of pulsars was chosen as the optimal fre- 
quency. To find the frequency that provided the largest 
number of pulsars, a cubic polynomial function was fitted 
through the points. The optimal centre frequencies for a 
pulsar survey inside and outside the Galactic plane are 
shown in Fig. [1] It can be seen that the optimal survey 
frequency decreases with lower sensitivity. The reason for 
this is that for lower sensitivity, the average distance to 
the detected pulsar is smaller. Thus, frequency dependent 
propagation effects are smaller, allowing lower frequen- 
cies to be used. It is surprising that the optimal centre 
frequencies for normal and millisecond pulsars are simi- 
lar, since scattering effects on the pulsar profile have a 
larger impact on the detection of millisecond pulsars than 
on the detection of normal pulsars. Thus it can be ex- 
pected that in the Galactic plane, higher observation fre- 
quencies would favour the detection of millisecond pulsars. 
However, milli second pulsars have a ste ep luminosity dis- 
tribution fe.g. ICordes fc Chernofi 1997 ). leading to many 
low-luminosity millisecond pulsars which are not detected 



4 



Smits et al.: Draft: Pulsar searches and timing with the SKA 



at high observation frequencies. It should be noted that 
with 500 MHz bandwidth, the lowest centre frequency is 
750 MHz for the 500-MHz dishes and 1.05 GHz for the 
800-MHz dishes. Once the optimal centre frequencies were 
determined, these restrictions were included in the simu- 
lations. We then performed four simulations correspond- 
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Fig. 1. Optimal frequency to survey the Galactic plane 
(|b| < 5°, |1| < 45°) or outside the Galactic plane (|b| > 
5°, |1| > 45°) as a function of sensitivity. The dotted hor- 
izontal lines show the lowest available centre frequency of 
the 500-MHz (bottom line) and 800-MHz (top line) dishes. 
Surveys with the AA were performed at a centre frequency 
of 650 MHz. 

ing to possible pulsar surveys. For the dishes they were 
performed at the optimal centre frequencies, taking the 
frequency limits into account. The frequency range of the 
AA was held fixed at 500 to 800 MHz. The simulations 
are: 

— An all-sky survey with the AA (including the Galactic 
plane) and a Galactic plane survey with the 800-MHz 
dishes. This shows the fraction of pulsars that can be 
found with configuration C. 

— An all-sky survey with just the AA. This shows the 
fraction of pulsars that will be found with configura- 
tion C, if we do not include the dishes in the survey. 

— An all-sky survey with the 500-MHz dishes. This shows 
the fraction of pulsars that can be found with config- 
urations A or B. 

— A Galactic plane survey with the 500-MHz dishes. This 
shows the fraction of pulsars that will be found with 
configuration A or B, without performing a survey out- 
side the Galactic plane. 

We assumed a bandwidth of 500 MHz for the dishes and 
300 MHz for the AA and an integration time of 1800 s. The 
simulations were limited to isolated pulsars only. A paper 
addressing the problems of the searching for and timing 
of binary pulsars is in preparation (Smits et al. in prep.) 

A Galactic centre survey with the SKA will allow 
for testing the conditions near a 3 x 10^ Mq black hole. 



Simulating such a survey requires a detailed modelling 
of t he scattering in the Galactic centre. Previous work 
fe.g.lLazio fc CordedlTiQl iBhat et al]l2004l : ICordes et al " 



20041 ) shows that frequencies above 10 GHz are required to 
penetrate the scattering screen surrounding the Galactic 
centre. These results can be used for a detailed study of 
a Galactic centre survey with the SKA. However, such a 
study is beyond the scope of this paper. Further results on 
this issue will be presented in J. Deneva et al. (in prep.) 

3.2.2. Simulation results 

Fig.lHshows the fraction of normal and millisecond pulsars 
that were detected in the simulation of 4 different surveys 
as a function of sensitivity. Only 1.5% of the pulsars could 
not be detected due to the declination of the SKA of -30°. 
As will be shown in Section[5l we expect that initially only 
the 1-km core of the SKA can be used for the pulsar sur- 
vey. This leads to a sensitivity of 0.1 SKA units, which 
is similar to having a fully steerable Arecibo-class tele- 
scope in the southern hemisphere. Performing an all-sky 
survey with only the 1-km core of the SKA, will detect 
about 14 000 normal pulsars out of a possible 30 000 and 
about 6 000 millisecond pulsars out of a possible 30 000, 
regardless of the SKA implementation. The total obser- 
vation time of the survey, however, does depend on the 
implementation. Using single-pixel feed dishes to perform 
an all-sky survey or a survey of just the Galactic plane 
will take 600 days and 30 days, respectively. With phased 
array feed dishes or an AA the total survey time can be 
much less, depending on the available computation power 
(see Section [6| . Because of the large natural FoV of the 
SKA and it's location in the southern hemisphere, the 1- 
km core of the SKA complements and extends work done 
by both Arecibo and the Five hundred meter Aperture 
Spherical Telescope (FAST), which is expected to be com- 
pleted around 2014. 



4. Timing 

All pulsars that will be found by the SKA need to be- 
come part of a regular timing program to find the inter- 
esting pulsars. Ideally, they need to be timed once every 
two weeks for about 6 months to a signal-to-noise ratio of 
about 9 (which is only a nominal target to characterise the 
new pulsars and to identify those sources which are inter- 
esting for high precision timingj^. The minimum timing 
duration is also restricted by the pulsar stabilisation time, 
which is the time it takes to obtain a stable integrated 
profile. Assuming that they will be timed using dishes, 
the total duration of timing all newly discovered pulsars 
just once depends on how the observations are performed. 



^ Usually, at least 12 months are needed to break the degen- 
eracies between positional uncertainties and pulsar spin-down. 
However, we assume that the imaging capabilities of the SKA 
will be used after the discovery to obtain an interferometric 
position, which speeds up the process. 
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Fig. 2. Fraction of normal (top) and millisecond (bottom) 
pulsars detected from different pulsar survey simulations 
as a function of sensitivity of the SKA. For both the 500- 
MHz and 800-MHz dishes, the surveys were performed at 
the optimal centre frequency, taking the frequency limits 
into account (see Fig.[T]). The frequency range of the AA 
was 500 to 800 MHz. 



i.e. at which points on the sky the telescope is pointed. 
To estimate the total duration of this timing process, we 
therefore suggest the following scheme that minimises the 
average observation time per pointing, by grouping dim 
pulsars together in one FoV. Let pfov be the angle of the 
FoV. 

1. All pulsars are placed on a long list, in order of in- 
creasing brightness. 

2. All pulsars within an angle of pfov from the first pulsar 
of the list (the dimmest pulsar) are placed on a short 
list. 

3. We now consider the dimmest pulsar on the short list. 
If possible, we restrict the pointing of the FOV to in- 
clude this pulsar and we remove this pulsar from both 
the short list and the long list. If the pointing of the 
FOV is already too restricted to include this pulsar, it 
is removed from the short list only. 



4. Step 3 is repeated until the short list is empty. The 
resulting pointing is stored. 

5. Step 2 through 4 are repeated until the long list is 
empty. 

This timing optimisation method saves about 30 to 50% 
of observation time over a simple grid approach. To com- 
pare timing performances, we used this scheme to esti- 
mate the time to obtain a single timing point of 14000 
pulsars from an all-sky survey to a signal-to-noise ratio of 
9 for each of the SKA-configurations, assuming usage of 
the full collecting area a nd a pulsar stabilisation time (see 
e.g. iHelfand et al.lll975h of 200 pulses. This leads to ob- 
servation times of 6 days for configuration A and 1.5 days 
for both configurations B and C. For configuration C the 
AA was used to time all the pulsars that were detected 
by the AA in the survey. These estimates do not include 
densely spaced 'resolving' observations that arc required 
to obtain the first phase coherent timing solution. Fig. [3] 
shows the simulated pulsar-sky and the pointings that re- 
sult from the timing optimisation method, assuming the 
20 deg^ FoV of configuration B. 




Fig. 3. Hammer projection of the Galaxy with 14000 de- 
tected pulsars from an SKA pulsar survey simulation. The 
horizontal axis corresponds to longitude and the vertical 
axis corresponds to latitude. The Galactic centre is located 
in the middle of the plot. The dashed lines running from 
top to bottom correspond to steps of 60° in longitude, 
whereas the dashed lines running from left to right corre- 
spond to steps of 30° in latitude. The circles indicate all 
pointings (700) that result from the timing optimisation 
method, assuming a FoV of 20deg^. 

To perform strong field tests of General Relativity and 
gravitational wave detection, as defined in the key science 
project, it is essential to time a large number of specific 
pulsars to signal-to-noise ratios of up to 100 on a regular 
basis. To estimate how much observation time is required 
for such high-precision timing, we have to distinguish be- 
tween timing millisecond pulsars (for the Pulsar Timing 
Array) and timing pulsars in a binary with a neutron star 
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or a black hole companion. We calculated the observation 
time for timing the 250 millisecond pulsars with the best 
signal-to-noise ratio's for the three configurations. We as- 
sumed that for timing purposes, the nearly full collecting 
area (out to several hundreds of kilometres from the core) 
of the SKA can be used. As a conservative estimate (sub- 
ject to current studies), we timed every pulsar for at least 
5 minutes to ensure a stable pulsar profile to minimise the 
error in the time of arrival. With the single-pixel feed 15- 
m dishes this takes about 20 hours, the 15-m dishes with 
phased array feeds take about 15 hours and the AA only 
takes about 6 hours. However, this assumes that the po- 
larisation purity after calibration of the AA is similar to 
that of the dishes, which might not be the case. To esti- 
mate the maximum required observation time for timing 
pulsars in a binary with a neutron star or black hole, we 
will assume a (somewhat optimistic) number of 200 such 
binaries that are potentially detectable in the Galaxy. We 
further assume that the characteristics of the pulsars in 
these binaries are similar to isolated pulsars. This leads 
to a detection of 90 compact binary pulsars. With the 
single-pixel feed 15-m dishes it takes about 2 days to time 
these binary pulsars, the 15-m dishes with phased array 
feeds take about 1.5 days and the A A only takes about 
8 hours (once again, assuming similar polarisation purity 
as the dishes). However, timing only the brightest 80% 
of the binary pulsars, would take about half the observa- 
tion time. It should be noted that these are only initial 
estimates to give an indication of the maximum required 
observation time for high-precision timing and to compare 
the performance of the different configurations. A paper 
containing a more detailed study on high-precision timing 
is in preparation (Smits et al. in prep.). 



5. Computational requirements 



A pulsar survey requires the coherent addition of the sig- 
nals from the elements in the core of the SKA and to form 
sufficient pencil beams to create the required FoV. These 
pencil beams will produce a large amount of data which 
will need to be analysed. We estimate the required com- 
putation power and data rates for 3 000 15-m dishes with 
a bandwidth of 500 MHz and 2 polarisations and for the 
AA with a frequency range from 500 to 800 MHz and 2 
polarisations. 



5.1. Beam forming 

Following ICordesI (|2007l l. we estimate the number of op- 
erations per second (ops), iVosb, to fill the entire FoV of a 
dish with pencil beams a^: 



Dr 



-Ddish 



(1) 



where Fc is the fraction of dishes inside the core, iVdish is 
the total number of dishes, A^poi is the number of polari- 
sations, B is the bandwidth, I?dish is the diameter of the 
dishes, -Dcorc is the diameter of the core and ^ is the factor 
by which the phased array feed extends the FoV. To fill 
the FoV of single-pixel 15-m dishes in a 1-km core requires 
4400 pencil beam. The number of operations per second 
to beam form a 1-km core containing 20% of 2 400 single- 
pixel 15-m dishes then becomes 2.2 x 10^^. The number of 
operations for beam forming the full FoV of 15-m dishes 
with a phased array feed are simply a factor ^ higher, 
corresponding to the increase in FoV. For the SKA im- 
plementation B, this corresponds to ^ = 31. The number 
of operations for beam forming can be reduced by beam 
forming in two stages, if we assume that the dishes in the 
core are positioned such that they can form sub-arrays 
that all have the same size. In the first stage, beam form- 
ing of the full FoV is performed for each sub-array. In the 
second stage the final pencil beams are formed by coher- 
ently adding the corresponding beams of each sub-array 
formed in the first stage. When the sub-arrays have the 
optimal size given by I?coro(^cA'dish)^^^*, as we show in 
the Appendix, the number of operations per second be- 



iVo: 



-Ddish 



(2) 



For the same parameters as above, the number of opera- 
tions per second becomes 2.0 x 10^^. For the beam forming 
of the AA, the calculations are slightly different. Initial 
beam forming will be performed at the stations them- 
selves, leading to beams equivalent to those of a 60-m dish 
for each station. Demanding a total FoV of 3 deg^ and as- 
suming a total collecting area of 500 000 m^ , the number 
of operations becomes: 



Nosh = F.NdishNpoiB ■ 3 D, 



2 

core max' 



(3) 



where c is the speed of light and A^dish ~ 177. For a 1- 
km core, this leads to a required computation power of 
1.4 X 10^** ops. Fig. m shows the number of operations per 
second for the beam forming as a function of core diameter 
for the 15-m dishes with single-pixel feed and the AA. 
To obtain the benchmark concentration as a continuous 



^ Note that for a phased array feed this function is actually 
frequency dependent, because the total FoV of the dish then 
becomes constant, yet the FoV of the pencil beams is given by 
(X/Dtei)'^. However, for the values of bandwidth and observa- 
tion frequency used in this paper, Eq.[T]is accurate. 
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Fig. 4. The number of operations per second required to 
perform the beam forming for the 15-m dishes and the 
AA. It is assumed that there are 2 400 15-m dishes and 
the bandwidth is 500 MHz. The FoV of 0.64 deg^ corre- 
sponds to fiUing the natural FoV of the dishes {£_ — 1) at 
a frequency of 1.4 GHz. For the AA we assume a FoV of 
3deg^ and a total coUecting area of 500 000 m^. The fre- 
quency range is 0.5 to 0.8 GHz. In all cases the number of 
polarisations is 2. The thick black line corresponds to the 
beam forming of the 15-m dishes in one stage. The thin 
black line corresponds to the beam forming of the 15-m 
dishes in two stages. The striped/dotted hue corresponds 
to the beam forming of the AA. 



function of core diameter between and 10 km, we used 
the following expression: 



a[l - exp(-6Dcoro)], 



(4) 



where a and b were tuned to the specifications of 20% and 
50% of the dishes within a 1-km and 5-km core, respec- 
tively, which leads to a = 0.56 and b — 0.45 x 10~^ m~^. 

As an alternative to beam forming by coherently 
adding the signals from the dishes up to a certain core 
size, it is also possible to perform the beam forming by 
incoherently adding the signals from sub-arrays. This pro- 
cess is similar to beam forming in two stages as mentioned 
above, except that in the second stage the beams that were 
formed in the first stage are added incoherently. This leads 
to a much larger beam size, which reduces the required 
computation power for beam forming significantly. It also 
reduces the total data rate and the required computation 
power for the data analysis. The drawback is that it re- 
duces the sensitivity of the telescope by a factor of ^/Nsa, 
where Ng^ is the number of sub-arrays, which will be sev- 
eral hundred. However, this can be partially compensated 
as this method allows utilising all elements that are placed 
in sub-arrays of the same size, which can possibly be as 
much as 5 times the collecting area of the core of the SKA. 



5.1.1. Data analysis 

There are two factors which impact the total data volume 
to be analysed. Firstly, the FoV will be split up in many 
pencil beams, each of which will need to be searched for 
pulsar signals. Secondly, the SKA will be able to see the 
majority of the sky, all of which we want to search for 
pulsars and pulsar binaries. There are two ways to achieve 
this. The first option is to analyse the data as it is received, 
immediately discarding the raw data after analysis. This 
requires the analysis to take place in real time. The second 
option is to store all the data from part of the survey and 
analyse them at any pace that we see fit before continuing 
with the next part of the survey. Both approaches pose 
serious technical challenges, which we will discuss here. 

First we consider the dishes for which we estimate the 
data rate of one pencil beam as 



V. 



dish 



1 B iVbits„ 

A^poi— ^ Bps, 



^samp 



8 



(5) 



where isamp is the sampling time, B is the bandwidth, 
Av is the frequency channel width, A'poi is the number of 
polarisations and A'bits is the number of bits used in the 
digitisation. Av can be estimated by demanding that the 
dispersion smearing within the frequency channel does not 
exceed the sampling time, which leads to: 



Ai/(GHz) = 



^samp(A^s)t/,^^i„(GHz 
8.3 X 103DM,nax 



(6) 



where i^min is the minimum (lowest) frequency in the ob- 
servation frequency band and DMmax is the maximum 
expected dispersion measure. For the dishes, the number 
of pencil beams to fill up the FoV can be estimated as 
C(£'core/£'dish)^- For the AA, the number of pencil beams 
becomes frequency dependent and for 3 deg^ FoV it is 
given by 



A^pbAA 



(7) 



However, a pencil beam needs to be centred on the same 
point on the sky over the frequency band to allow for sum- 
ming over frequency. Thus, even though the pencil beams 
become larger towards lower frequency, the number of pen- 
cil beams needs to remain constant and is determined by 
the highest frequency. The total data rate from the AA 
can then be estimated as 



■Daa = 



1 



^samp 



B iVbits„ 
Apoi— 7— 3 



vi80cy 



D. 



core max 



Bps. (8) 



Fig. El shows the data rate from a pulsar survey for the 
15-m dishes and the A A as a function of core diame- 
ter, assuming tgamp^lOO^s, DMmax^lOOOcm^"^ pc for the 
dishes and DMinax^SOOcm^"^ pc for the AA, 7Vpoi=l (sum 
of 2 polarisations) and 2 bit digitisation. The AA was as- 
sumed to operate on the full frequency range of 0.5 to 
0.8 GHz. The frequency range for the dishes was assumed 
to be 1 to 1.5 GHz. For a 1-km core, the data rate for the 
15-m dishes with single-pixel feed is 4.7 x 10^^ bytes per 
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Fig. 5. Data rate from pulsar surveys using the 15-m 
dishes or the AA, as a function of core diameter. 

second and for the AA it is 1.6 x 10^^ bytes per second. 
Both values scale linearly with the FoV. We can estimate 
the number of operations to search one 'pencil beam' for 
accelerated periodic sources as one Fourier Transform of 
all the samples in the observation, for each trial DM-value 
and for each trial acceleration: 



^oa = A^DMA^acc X 57Vsamp logs ( A^'samp) , 



(9) 



Adm — 



4150DM,„ax(^^-?n(GHz) 



7-2 

max 



(GHz)) 



t. 



samp 



(ms) 



(10) 



For a 1-km core, the required computation power be- 
comes 1.2x 10^^ and 4. Ox 10^^ ops for the 15-m dishes and 
the AA, respectively. These values are lower estimates of 
the actual number of operations, as there will be signif- 
icant contributions from dedispersion, harmonic folding 
and possibly other processes that will contribute to the 
total number of operations. These additional processes are 
all on the order of Agamp and we expect their total con- 
tribution to be similar to one Fourier Transform at most. 

Alternatively, we can consider storing the data from a 
part of the survey and analyse it at a much slower rate 
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Fig. 6. Operations per second required to perform a real 
time analysis of pulsar surveys using the 15-m dishes or 
the A A, as a function of core diameter. 



where Adm is the number of trial DM- values, Aacc is the S 
number of trial accelerations which scales with the square § 
of Asamp, which is the number of samples in one observa- s 
tion. Thus, the number of operations for the analysis scales 
as Ag^amp log(^samp)- This mcaus that increasing the ob- 
servation length is computationally very expensive. Once 
again, for the dishes the number of pencil beams to fill up 
the FoV can be estimated as ^(-Dcorc/-Ddish)^ and for the 
AA the number of pencil beams is given by Eq. [3 with 
'^^'^max^O.S GHz. Fig. [6] shows the number of operations 
per second required to perform a real time analysis of a 
pulsar survey with the 15-m dishes and the AA as a func- 
tion of core diameter, assuming 100 trial accelerations, a 
sampling time of 100 /iS and an observation time of 1 800 s. 
Adm was taken to be equal to the time shift due to the 
maximum dispersion divided by the sampling time, which 
is equal to 
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Fig. 7. Total amount of data from an all-sky survey 
(35 000 deg2) and a survey of the Galactic plane (900 deg^) 
as a function of core diameter, assuming a frequency range 
of 1 to 1.5 GHz. A 1-km core would correspond to 20% of 
the collecting area of the SKA. 



(see further discussion in Section 16. 2p . From Eq. [5] we can 
estimate the total amount of data from a survey. Fig. [7] 
shows the total amount of data from an all-sky survey 
and a survey of the Galactic plane as a function of core 
diameter, assuming a frequency range of 1 to 1.5 GHz, 
<samp=100/is, DMmax=1000cm-3pc, Apoi=l, an obser- 
vation time of 1 800 s and 2 bits digitisation. Assuming a 
storage capacity of 1 exa-byte, an all sky survey would 
have to be split up in about 40 parts. The total time of 
the survey would then depend on the speed at which the 
analysis takes place. 
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Center Frequency (GHz) 

Fig. 8. The frequency dependence of the figure of merit. 
This function expresses the speed of the survey for a fixed 
sensitivity of the telescope. 



5.1.2. Survey speed 

To express the trade-off between the SKA design and sur- 
vey speed, we define the figure of merit for dishes as (see 
appendix for derivation) 



^dish-Ddish 



,2q-2 



(11) 



sky V 400 MHz 



where rj is the aperture efficiency, ^ is the factor by which 
the phased array feed extends the FoV, B is the band- 
width, A'dish is the number of dishes used in the survey, 
Z?dish is the diameter of the dish, Ttci is the system tem- 
perature of the telescope (K), T^^^ is the sky temperature 
(K) at 400 MHz (Tgky scales as v is the observation 

frequency, and a is the av erage spectral index of radio pul- 



sars, which is about -1.6 (jLorimer et al.lll995f ). The figure 



of merit shown here expresses the speed of the survey for 
a required sensitivity of the telescope and is identical to 
that for a continuum survey of steady sources. It does not 
include the loss of sensitivity associated with pulse broad- 
ening effects and does not express the impact of dedisper- 
sion and periodicity search on the survey time. 

Fig. [5] shows the strong dependence of the figure of 
merit on the observation frequency. The slope at high 
frequencies is due to the spectral index of pulsars and 
the decrease of the FoV. At low frequency the figure of 
merit becomes flat as the sky temperature begins to dom- 
inate the system temperature. In practise, pulsars become 
harder to detect at low frequencies where scattering ef- 
fects broaden the pulsar profile. Since the scattering de- 
pends on the location on the sky, it was not included 
in the merit function. Thus, at low frequencies the fig- 
ure of merit is overestimated. However, scattering effects 
are only significant for pulsars located in the Galactic 
plane. For non-mill isecond pulsars, low frequenci es are 



very favourable (see Ivan Leeuwen fc: Stappersll2008( ) 



The total observation time of a survey is given by: 



total 



-t, 



FoV 



point 



(12) 



where flgur is the total solid angle of sky covered by the 
survey, fipoV is the total FoV of the dishes and ipoint is 
the observation time per pointing. 

6. Suggestions 

A n estimate of the computational power by 2015 is given 
bv lCornwelll (|2005f) to be lOPflop for $100 M. To estimate 
what the SKA can achieve for pulsar searches and tim- 
ing, we will assume that several 10^^ ops are available for 
beam forming and data analysis when this survey will be 
performed, which will be after 2020. 

We consider three scenarios, two with full coherent 
beam forming, where the analysis takes place either in 
real time or off-line and one scenario where part of the 
beam forming takes place incoherently. 

6.1. Real time analysis 

First we assume that the analysis takes place real time. 
The immediate benefits are that no time is lost for the data 
analysis and that only reduced data needs to be stored. 
One drawback, however, is that the data can only be pro- 
cessed once. Experience with current pulsar surveys show 
that it is extremely advantageous to have multiple passes 



Faulkner et al 



at pu lsar survey data analysis (see e.g. 
20051) . 

As calculated in Section I5.1.1[ for a 1-km core the 
data rate from the dishes with a single-pixel feed is 
4.7 X lO^^bytes/s. The required computational power for 
real time analysis is 1.2 x 10^^ ops. Both values scale lin- 
early with FoV. Only with SKA implementation B, where 
the dishes have phased array feeds, can the FoV of the 
dishes be increased beyond 0.64 deg^. Using the 1-km core 
of the AA and using a FoV of 3 deg^ leads to a data rate 
of 1.6 X 10^^ bytes/s and requires a computation power of 
4.0 X 10^^ ops to analyse the data in real time. 

SKA implementation A suggests the all-sky survey to 
be performed with single-pixel 15-m dishes. At low fre- 
quencies the FoV of the dishes is about 1 deg^ . The entire 
survey will then take about 600 days. (This value assumes 
100% telescope time. In practice this survey might take 
up to 5 years to perform). In implementation B a larger 
FoV can be used. The survey time of 600 days then scales 
inversely with the FoV, however the data rates and re- 
quired computation power go up linearly with the FoV. 
Implementation C has the benefit of the AA for which an 
all-sky survey takes only about 200 days and the survey 
of the Galactic plane with single-pixel 15-m dishes takes 
about 30 days. The data rates and required computation 
power for the AA are larger than for the dishes. When 
we take this into account the total survey time for all 
implementations become similar, except that implemen- 
tations B and C allow for a faster survey and the AA in 
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implementation C might be use d for other observations 
simultaneously (see ICorde£ll2007i for using the SKA as a 
synoptic survey telescope). 

6.2. Off-line analysis 

The off-line analysis requires the full data from an ob- 
servation to be stored. Thus, the data rates from Fig. [5] 
become the rates at which data is written to a storage de- 
vice. When the maximum data storage has been reached, 
the observation can be stopped and the data analysis can 
be performed at any pace suitable. After the data analysis 
has been completed, the next observation can be run. This 
allows a trade-off between computation cost and survey 
speed. A drawback is that, depending on available future 
storage solutions, the total survey time may easily become 
more than a decade for any SKA implementation. 

6.3. Incoherent beam forming 

For real time analysis, restrictions on the computational 
power and data rate will likely limit the FoV to about 
3 deg^. The survey time will then become at least 200 days. 
Off-line analysis requires the huge data rates to be written 
to a storage device and will increase the survey time by 
a significant factor. We therefore consider the possibility 
of incoherent addition of the sub arrays, as described in 
Section 15.11 

We will assume that all the elements in the SKA are 
placed such that their signals can be combined to form 
sub-arrays (or stations in the case of the AA) of 60 me tres 
in diameter (similar to Ivan Leeuwen fc Stappersll2008L for 
LOFAR) . The pencil beams that are formed by coherently 
adding the signals from the sub-array elements are about 
17 times larger than the pencil beams from coherently 
beam forming the 1-km core. Thus the number of pencil 
beams required to fill the FoV becomes about 280 times 
smaller. The computational requirements for beam form- 
ing, data analysis as well as the data rates go down linearly 
by this factor. This method of beam forming utilises the 
full collecting area of the SKA, which increases the sensi- 
tivity by a factor of 5. However, because the sub-arrays are 
added incoherently, the sensitivity decreases by the square 
root of the number of sub-arrays. For the AA there will 
be about 180 stations. Thus the sensitivity goes down by 
a factor of -\/180/5 « 2.7. Assuming that the dishes are 
placed in sub-arrays with a 20% filling factor, the num- 
ber of sub-arrays will be between 600 and 900, depending 
on the SKA implementation. For incoherent beam forming 
the sensitivity will then go down by a factor between 5 and 

6. Applying this to an all-sky survey leads to a detection 
of about 20% of all the detectable pulsars. 

7. Summary 

We have investigated the pulsar yields for different SKA 
configurations by simulating SKA pulsar surveys for dif- 
ferent collecting areas and centre frequencies. For the 



Galactic plane, the optimal centre frequency lies just 
above 1 GHz. Outside the Galactic plane the optimal cen- 
tre frequency lies between 600 and 900 MHz, depending on 
the collecting area. Combining the dishes and the AA to 
perform a survey inside the Galactic plane and outside the 
Galactic plane, respectively, would result in detection of 
roughly 14000 normal pulsars and 6 000 millisecond pul- 
sars. As Fig.[2]shows, increasing the collecting area within 
the 1-km core would improve the detected number of pul- 
sars significantly. The number of detected pulsars scales 
roughly as where Fc is the fraction of collecting area 
in the core. 

We also describe a simple strategy for follow-up timing 
observations and find that, depending on the configura- 
tion, it would take 1-6 days to obtain a single timing point 
for 14 000 normal pulsars. Configuration C, containing the 
AA, will be of great benefit here. Obtaining a single timing 
point for the high-precision timing projects of the SKA, 
will take less than 14 hours, 2 days, or 3 days, depending 
on the configuration. Once again, the presence of the AA 
allows for the fastest timing. In the case of high-precision 
timing, however, this assumes that the polarisation purity 
of the AA is similar to that of the dishes. 

Performing a pulsar survey with the SKA requires the 
coherent addition of the signals of the individual elements, 
forming sufficient pencil beams as to fill the entire FoV 
of a single dish. Because of the extreme computational 
requirements that arise due to large baselines, it is not 
possible to combine the signals of all of the elements in 
the SKA. Rather, only the core of the SKA can be used 
for a pulsar survey. We have derived the computational 
requirements to perform such beam forming as a function 
of core-diameter for the 15-m dishes and the AA. For a 1- 
km core the requirements are 2.2x10^^ and 1.4x10^'' ops 
for the 15-m dishes and the AA, respectively. When the 
dishes are placed such that they can form sub-arrays, the 
computational requirement for beam forming goes down 
significantly. We have also calculated the data rates and 
the computational requirements for applying a search- 
algorithm to the data to find binary pulsars. Both limit 
the usage of the SKA for pulsar searches to a core of about 
1 km. 
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Appendix 

Beam forming of the dishes 

Cordes (2007) estimated the number of operations to perform 
the beam forming for a pulsar survey as one operation (corre- 
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spending to the phase shift of one element) for each dish and 
each polarisation at the Nyquist frequency. This needs to be 
performed for each of the pencil beams that fill the total FoV. 
This leads to the total number of operations per second 



A^osb = 



FciVdishiVpoiB f ^) , 



(13) 



where Fc is the fraction of dishes inside the core, A'dish is the 
total number of dishes, A^poi is the number of polarisations, 
B is the bandwidth, Ddish is the diameter of the dishes and 
-Dcoro is the diameter of the core. The fraction (Dcorc/Odish)^ 
is equal to the required number of pencil beams to fill up the 
FoV. If we assume that the dishes in the core are positioned 
such that they can form sub-arrays, beam forming can take 
place in two stages. In the first stage beam forming of the full 
FoV is performed for each sub-array. In the second stage the 
final pencil beams are formed by coherently adding the cor- 
responding beams of each sub-array formed in the first stage. 
This leads to the following expression for the number of oper- 
ations required for beam forming in two stages: 



/ -Dsa \^ 

A'osb2 ~ NsaNdishsaNpolB I -— J 

V i^dish / 

i^j 



+ 



(14) 



where Nsa is the number of sub-arrays in the core, A'dishsa is the 
number of dishes in one sub-array and Dsn is the diameter of 
the sub-arrays. Substituting iVdishsa = -F'cA''dish/A''sa and A''sa = 
[Dcorc/Dsa)^ yields 



V i^dish / 



(16) 



which has a minimum at Dsn = £'corc(^cA'dish)~^''^. At this 
minimum, the number of operations for beam forming in 2 
stages is: 



NZ^2 = 2VFeiVdishArpolB 



/-Dcoro y 
V -Ddish / 



(16) 



Beam forming of the AA 

For the AA, initial beam forming will be performed at each 
station; this will lead to beams equivalent to those of a 60- 
m dish. Assuming a collecting area of 500 000 m^, this leads 
to about 177 stations. The FoV of 1 pencil beam is given 
by (A/Dcoro)^ = (c/i/Dcore)^ steradians, where c is the speed 
of light. Thus, to cover 3deg^ the required number of pencil 
beams is (37r/180)^/(c/i^-Dcoro)^. Within the frequency band, 
the number of pencil beams needs to remain identical as to en- 
able summing over frequency. In order to have a full coverage 
of the 3 deg^ , the required number of pencil beams is therefore 
obtained by setting u to i/max- The number of operations per 
second to perform the beam forming of the AA is then given 
by (see also Eq. [3} 



A^osb = FciVdishA^poiB • 3 



180c 



core max • 



(17) 



Derivation of figure of merit 

Neglecting numerical factors and physical constants, the figure 
of merit expresses the speed of a survey for a given sensitivity 
of the telescope. It therefore scales with aperture efficiency (rf), 
the FoV, the bandwidth (B), the square of the collecting area 
( A) and the reciprocal of the square of the system temperature 
(Tsys). For a pulsar survey it also scales with the square of 
the flux of pulsars, which scales as v"' , with v the observation 
frequency and a the average spectral index of radio pulsars, 
which is about -1.6. Thus: 



M = 7?BFoV 



Tsys 



(18) 



The FoV is given by £,{^0^^^^^^ , where ^ is the factor by which 
the phased array feed extends the FoV and Ddish is the di- 
ameter of the antennas. This is equal to ^icu'^ D~^.^^^^^ . The 
collecting area can be expressed as A = A^'dishTrD^igj-j/l, where 
A'^dish is the number of dishes used in the survey. Tsys is the 
sum of the telescope temperature and the sky temperature, 
the latter scaling with frequency as v~'^'^ . Putting all this in 
Eg. 1181 and removing all constants, yields: 



A^'dish-Ddish 



sky V400 



a V 

MHz ) 



(19) 



M = -qS^B 



Glossary 

The following list summarises the meaning of the most impor- 
tant parameters in this paper. 

Tsys: system temperature (including contributions from re- 
ceiver and sky background. 
r4°°: sky temperature at 400 MHz. 

FoV: Field of view, size of the sky area instantaneously covered 

by the telescope. 
T): aperture efficiency for dishes. 
Ddish: diameter of a single dish. 

Doi: distance between the furthest elements that are used in 
an observation. For any pulsar survey in this paper this 
parameter is equal to Dcore- 

Dcoro: diameter of the core region. 

Aosb: number of operations per second required for beam 

forming. 
Adish: total number of dishes. 

Fc'. fraction of the number of dishes inside the core. 
Apoi: number of polarisations. 
B: observing bandwidth. 

5: factor by which the phased array feed expands the FoV. 
Note that for a phased array feed this factor is frequency 
dependent. 

Vmax- maximum observing frequency. 

a, b: parameters appearing in Eq. [l] 

©dish: survey data rate for one pencil beam. 

T>AA'- total survey data rate for AAs. 

isamp: sampling time for pulsar survey. 

Av: channel bandwidth for pulsar survey. 

Abits: number of bits used for digitisation of survey data. 

Vmin- lowest frequency used for pulsar survey. 

DMmax: maximum dispersion measure to be searched or ex- 
pected. 

Adm: number of trial DM- values. 

A^s, 



number of samples in each survey time series. 
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Nacc'- number of trial accelerations to be searched. 
-^pbAA: number of pencil beams for A As. 

A^oa: number of operations to search one pencil beam in accel- 
eration searches. 
M: survey-speed figure-of-merit for searches with dishes, 
ttotai: total observing time for a survey, 
tpoint: observing time per survey pointing. 
Qaur'- total sky (solid angle) for survey. 
Qfov: total FoV for dishes. 
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